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Torsion-rotation transitions in molecules exhibiting hindered internal rotation possess enhanced 
sensitivities to a variation of the proton-to-electron mass ratio. This enhancement occurs due to 
a cancellation of energies associated with the torsional and rotational degrees of freedom of the 
molecule. This effect occurs generally in every internal rotor molecule, but is exceptionally large 
in methanol. In this paper we calculate the sensitivity coefficients of methyl mercaptan, the thiol 
analogue of methanol. The obtained sensitivity coefficients in this molecule range from K^j, = —14.8 
to +12.2 for transitions with a lower-level excitation energy below 10cm~^. 
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I. INTRODUCTION 

Physical theories extending the Standard Model of 
particle physics have presented scenarios that allow for 
spatial-temporal variations of the constants of nature [1] . 
Since the initial findings of a possible variation of the fine 
structure constant by Webb et al. [2] there has arisen a 
great activity in search for signatures of such variations. 
Studies aimed at detecting a possible drift of a fundamen- 
tal constant on a cosmological time scale focus mainly on 
the fine structure constant, a, [3-6] and the proton-to- 
electron mass ratio, /z [7-9]. A variation of a or /i will 
manifest itself as a change in the spectrum of atoms and 
molecules, since such a variation may induce a shift in the 
position of a spectral line. Not all lines will shift in the 
same amount or direction. The response of a transition 
to a variation oi a ov fi is characterized by its sensitivity 
coefficient, K^ or K^, respectively, which is defined as 
the proportionality constant between the fractional fre- 
quency shift of the transition, lS.v /v and the fractional 
shift in a or /i. 



= Kx —rr- , 

V X 



with X = 



a, ^. 



(1) 



The search for /i variation on a cosmological time scale 
has been made operational by comparing optical transi- 
tions of molecular hydrogen (Hj) in high-redshifted ob- 
jects with accurate laboratory measurements [7]. These 
investigations have yielded a limit at the level of A/z//x < 
10~^ for look-back times of 12 billion years [8, 9]. The 
transitions in Hj that were used to obtain this result 
possess sensitivity coefficients that range from —0.05 to 
-1-0.02. Inversion transitions of ammonia (NHg) were 
found to be ^100 times more sensitive to /i-variation 
than Hg transitions [10, 11]. Astronomical observations 
of NHg, in the microwave or radio range of the electro- 



magnetic spectrum, led to stringent \a constraints at the 
level of (I.0± 4.7) xlO~^ [12] and (3.5± 1.2) x 10"^ [13]. 
Soon thereafter it was realized that the large number 
of degrees of freedom that exist in even the simplest 
polyatomic molecules can result in large enhancements 
of the sensitivity coefficients for a possible drift in /x. 
These enhancements occur for transitions between near- 
degenerate levels that each have a different dependence 
on /i. For instance, it was found that mixed inversion- 
rotation transitions in HgO"*" have sensitivity coefficients 
ranging from K^ = —9 to -1-5.7 [14], while the Renner- 
Teller interaction in /-CHg results in sensitivity coeffi- 
cients ranging from K^j^ = —53 to -1-742 [15]. Mixed 
torsion-wagging-rotation transitions in methylamine dis- 
play sensitivity coefficients ranging from K^ = —19.1 to 
-0.75 [16]. 

In the context of astrophysical searches, methanol [17, 
18] is the target species of choice since it possesses sen- 
sitive transitions at low excitation energy and has been 
observed at high redshift [19]. In a recent study, Bagdon- 
aite et al. used four transitions in methanol to constrain 
A/x//x at (0.0 ± 1.0) X 10"'^ at a look-back time of 7 bil- 
lion years [20]. Methyl mercaptan (CHgSH) is the sulfur 
analogue of methanol and might therefore posses transi- 
tions that have large sensitivity coefficients to a variation 
of jjL [21]. Although methyl mercaptan has thus far only 
been detected in our local galaxy [22, 23], recent advances 
in radio telescopes have greatly increased the number of 
detected molecular species at high redshift. It is therefore 
relevant to have a list available that contains the sensi- 
tivity coefficients of transitions in methyl mercaptan that 
might be observed in the interstellar medium. 

The recent terahertz and far-infrared study of the nor- 
mal isotopologue (^^CHg'^SH) of methyl mercaptan by 
Xu et al. [24] has resulted in a complete list of the 
molecular parameters for this molecule. In this paper 
we use the results of Xu et al. [24] and the scaling rela- 



tions of the molecular parameters that were derived in 
Jansen et al. [17, 21] to calculate the sensitivity coeffi- 
cients of methyl mercaptan. 



II. STRUCTURE OF METHYL MERCAPTAN 

Methyl mercaptan or methanethiol, depicted on the 
right-hand side of Fig. 1, consists of a thiol group (SH) 
attached to a methyl group (CHg) and is thus the sulfur 
analogue of methanol. The CS sigma bond connecting 
the two parts of the molecule is flexible, allowing the 
methyl group to rotate with respect to the thiol group. 
As in the case of methanol, this rotation is not free 
but hindered by a threefold potential barrier with min- 
ima and maxima that correspond to the staggered and 
eclipsed configuration of the molecule, respectively. For 
the lowest energy levels this relative or internal rotation 
is classically forbidden and only occurs due to quantum 
mechanical tunneling of the hydrogen atoms. As a con- 
sequence of this tunneling each rotational level splits into 
three levels that are labeled according to their symme- 
try as A or E^ as can be seen in on the left-hand side of 
Fig. 1. 

The lowest energy levels of CHgSH are shown in the 
left and right panel of Fig. 2. The A and E species can be 
considered as two different molecular species in the same 
sense as ortho- and para-ammonia, respectively. The ar- 
rangement of energy levels within a symmetry state re- 
sembles that of a prolate symmetric top, with the differ- 
ence being that every K ladder obtains a small energy 
offset due to the K dependent tunneling splitting. 

As a consequence, certain states in neighboring K lad- 
ders may become near degenerate which results in a large 
enhancement of the sensitivity coefficients K^^ for transi- 
tions connecting these states [17, 18]. 



III. SENSITIVITY COEFFICIENTS 
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Figure 1. (color online) Variation of the potential energy 
of methyl mercaptan with the relative rotation 7 of the SH- 
group with respect to the methyl group about the molecular 
axis. Shown are the J = 1, jJ^j = 1 energies of the lowest 
torsion-vibrational levels. The splitting between the different 
symmetry levels is due to tunneling through the potential 
barriers. The j4-symmetry species are split further due to the 
asymmetry of the molecule (-ft'-splitting). 



with /i. These scaling relations were obtained in a pre- 
vious study [17], while a more detailed discussion on the 
derivation can be found in Ref. [21]. 

Table I lists calculated transition frequencies and sen- 
sitivity coefficients in methyl mercaptan with an excita- 
tion energy less than 10cm~^, that is, those transitions 
most relevant for astrophysical searches. It can be seen 
that several of these transitions display relatively large 
sensitivity coefficients. 

In the last column of Table I the results of the approx- 
imate model of Ref. [21] are listed. This "toy" model is 
derived for molecules that exhibit hindered internal ro- 
tation and contain a C^^ symmetry group. The model 
decomposes the energy of the molecule into a pure rota- 
tional and a pure torsional part. The rotational part is 
approximated by the well-known expression for the rota- 
tional energy levels of a slightly asymmetric top 



The energy levels of methyl mercaptan have been cal- 
culated using an adapted version of the belgi code [25] . 
This code was modified and improved by Xu et al. [26] 
in a number of ways useful for treating the large datasets 
available for the methanol molecule, but the general ap- 
proach has not been changed compared to the initial 
code. The present version of the code has been sped 
up compared to the original version, and also a substan- 
tial number of higher order parameters has been added. 
Using the set of 79 molecular parameters of methyl mer- 
captan obtained by Xu et al., the lower energy levels are 
found with an accuracy < 100 kHz. 

In order to calculate the K^^ coefficients for the differ- 
ent transitions in methyl mercaptan, the energy of each 
level and its dependence on fi has to be obtained. This 
translates into knowing the values of the molecular con- 
stants that go into belgi and how these constants scale 



ErotiJ, K) = -iB + C)J{J + l)+{A 



B + C 



K' 



(2) 
C = 
and c 



with A = 3.428 cm- 1, B = 0.432 cm^^, and 
0.413 cm~^ the rotational constants along the a, b 
axis of the molecule, respectively. The torsional energy 
contribution is approximated by a Fourier expansion as 



Eto.s{K) = F 



Qq + fli cos 



27r 



[pK 



(3) 



with F = 15.040 cm^^ the constant of the internal ro- 
tation, p — 0.652 a dimensionless constant reflecting 
the coupling between internal and overall rotation, and 
cr = 0,±1 a constant relating to the torsional symmetry. 
The expansion coefficients aq and ai depend on the shape 
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Figure 2. Energy of the lowest rotational levels in the torsion-vibrational ground state {ut = 0) of methyl mercaptan 
( CHg'^SH) [24]. The levels are denoted by Jk (indicated on the left side of each level). For the A levels the so-called 
parity quantum number (+/— ) is also used. The panel on the left displays the A state levels, whereas the panel on the right 
displays the E state levels. High sensitivities are expected for transitions that connect near degenerate levels with different K. 




441.442 cm^^ the height of the barrier, is a measure of 
the effective potential. The torsional energy for methyl 
mercaptan is plotted in Fig. 3 as a function of K. Note 
that the torsional splitting between the A and E levels in 
the K — state of methyl mercaptan is 0.99 cm~^, and 
thus an order of magnitude smaller than in methanol, 
which has a torsional splitting of 9.1 cm^^. As a conse- 
quence, the amount of energy that can be cancelled in 
methyl mercaptan will be less than in methanol. 

Finally, the sensitivity coefficient of the transition is 
obtained from 



Figure 3. (color online) Torsional energies in the ground tor- 
sional state {ft = 0) of methyl mercaptan obtained with belgi 
for A (solid circles), El (open circles), and E2 (open trian- 
gles) levels as function of K. The solid, dashed, and dotted 
curves are fits to Eq. (3) for A, El, and E2 states. Note that 
only integer values of K have physical meaning. 



of the torsional potential. Since we are mainly interested 
in the torsional energy difference ao cancels, ai follows 
from 



ai = Ais 

with Ai = -5.296, Bi = lAU, 
The dimensionless parameter s 



-Ci^s 



(4) 



and Ci = 2.120 [21]. 
= W^/%F, with V-i = 



Jftoy 



^rot 



A-B,, 



K^^'^^Eu 



AE„ 



AE,. 



(5) 



Note that rather than using A-Brot + AE'tors, we chose to 
use the experimental energy difference between the lev- 
els, hf, in order to account for the slight asymmetry of 
the molecule. The sensitivity of a pure rotational transi- 
tion is -K^™* = — 1, whereas the sensitivity of a pure tor- 
sional transition is given by Kj^'-'^ — {Bi — l) — ^Ciy/s — 
—3.7 [21]. The sensitivity coefficients of this simple model 
are seen to agree well with the results obtained by a di- 
agonalization of the full molecular Hamiltonian, which 
reflects the robustness of the obtained results. 



Table I. Transition frequencies and sensitivity coefficients in metliyl mercaptan with a lower-level excitation energy of less than 
fOcm"^, calculated with belgi using the molecular constants of Ref. [24] (fifth column) and the toy model of Ref. [21] (sixth 
column). The third column lists the transition strength multiplied by the electric dipole moment, /le, squared, while the fourth 
column lists the excitation energy of the lower level in Kelvin. 
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IV. CONCLUSIONS 

In this paper we have calculated sensitivity coefficients 
for transitions between low lying rotation levels in methyl 
mercaptan. The reported sensitivities span a range from 
K^ = —14.8 to +12.2 and can therefore be used to search 
for variation of fx in methyl mercaptan only. Although, 
thus far methyl mercaptan has only been detected in our 
local galaxy, it is our hope that the advanced spectral 
coverage, resolution, and sensitivity of the new genera- 
tion of radio telescopes such as ALMA (Atacama Large 
Millimeter Array) will result in the detection of this 
molecule at high redshift. Note that the comprehensive 



line list of accurate rest frequencies for methyl mercap- 
tan obtained by Xu et al. [24] should also alleviate this 
search. 
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